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We construct an exact dynamical mean field theory for nonresonant inelastic light scattering in 
the infinite-dimensional Falicov-Kimball model, which can be tuned through a quantum critical 
metal-insulator transition. Due to the projection of the polarization orientations onto different 
regions of the Brillouin zone and due to the transfer of energy and momentum from the light to 
the strongly correlated charge excitations, the nature of the dynamics can be naturally interpreted 
as strongly temperature-dependent low-energy particle-hole excitations and weakly temperature- 
dependent high-energy charge transfer excitations which depend delicately on the electronic corre- 
lations. These results can be used to give important information concerning the evolution of charge 
dynamics in different regions of the Brillouin zone. 

PACS numbers: 78.70.Ck, 72.80.Sk, 78.66.Nk, 71.30.-fh, 71.27.-fa 



I. INTRODUCTION 

Inelastic x-ray scattering (with meV accuracy over a 
spectral range of several eVs) has improved significantly 
over the past few years due to the increased photon flux 
of third generation synchrotron gourcesiiSi^iiiSiSiLSiSiiSiii. 
The large cross-section of light-coupled probes (as com- 
pared to neutron scattering, for example) allows for a 
systematic study of the dispersive charge dynamics in a 
wide dynamical range (q, f2) in solids and fluids. It has 
opened an additional window to study correlation effects 
on phonons^, plasmons^, quasiparticles'', charge transfer 
excitation3^iSiL2i2ii£, and orbital excitationsii. One par- 
ticular point of interest has been the study of the evo- 
lution of strongly correlated systems as some parameter 
of the system, such as the electron density, is varied by 
doping or pressure. While many single-particle proper- 
ties have been studied via angle-resolved photoemission, 
important questions concerning the evolution of the un- 
occupied states are now directly accessible via inelastic 
x-ray scattering. 

Recent experiments have focused on a number of 
correlated (Mott) insulators such as La2Cu04 and 
Sr2Cu02Cl2^ Ca2Cu02Cl2-7, NaVzOgS, NdsCuOA and 
the one-dimensional insulators Sr2Cu03 and Sr2Cu02iS. 
The measurements have revealed dispersive high-energy 
and low-energy excitations which have been identified 
with a photon-induced charge transfer between different 
atomic orbitals, or with transitions from the lower to 
the upper Hubbard band across an effective q-dependent 
Mott gap. 

More recent measurements have begun to appear 
in materials doped from their parent Mott insulating 
phasesi^. However, the theoretical development of in- 
elastic x-ray scattering in strongly correlated metals and 
insulators is just starting to formiiiii^iiSiii. Of partic- 
ular interest is a determination of how the upper and 
lower Hubbard bands and consequently the Mott gap 



evolve with correlations. As experiments reach greater 
and greater resolution, it will shortly be possible to track 
the evolution of electronic correlations from strongly cor- 
related insulators to strongly and then weakly correlated 
metals. The purpose of this contribution is to investi- 
gate such a theory for inelastic x-ray scattering. In par- 
ticular, we develop an exact dynamical mean field the- 
ory for nonresonant inelastic light scattering in a system 
which can be tuned across the quantum critical point of 
a metal-insulator transition. We calculate the inelastic 
x-ray cross-section on both sides of the transition and 
near the critical point. 

The outline of this manuscript is as follows: in Section 
II we develop the general formalism for nonresonant in- 
elastic x-ray scattering and review simple physical ideas 
for weakly correlated metals. In Section HI we present 
the specific formalism for calculating the x-ray response 
in the Falicov-Kimball model in the limit of large spatial 
dimensions, and in Section IV we present the numerical 
results. Lastly, we summarized our results and discuss 
them in light of recent measurements in Section V. This 
paper expands on results for the insulating phassi^ to 
consider metals and materials close to the metal-insulator 
transition. 



II. FORMALISM 

A. Non-resonant response 

Light can scatter off of many different excitations in a 
system, but here we focus on the inelastic scattering of x- 
rays from electrons. X-rays, unlike optical photons, can 
exchange both energy and momentum when they scat- 
ter with a solid. The scattering occurs as light creates 
charge fluctuations in different locations of the Brillouin 
zone (BZ). These charge fluctuations are classified as ei- 
ther isotropic charge fluctuations or anisotropic charge 
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fluctuations (which vanish when averaging their k-space 
variation over the BZ). The way in which the charge fluc- 
tuations are created is dictated by the polarization ori- 
entation of the incoming and outgoing photons set by 
the scattering geometry. These polarization orientations 
transform according to the operations of the point-group 
symmetry of the crystal, and therefore so must the charge 
fluctuations that they create. It is through this mecha- 
nism that the charge excitations in different regions of 
the BZ can be systematically selected and explored via 
inelastic light scattering. 

These charge fluctuations relax by internal scattering 
processes, such as due to impurities or Coulomb scatter- 
ing, and finally via the re-emission of photons; inelastic 
x-ray scattering probes these relaxation processes at dif- 
ferent regions of the BZ and at different transferred en- 
ergies. An important distinction between isotropic and 
anisotropic charge fluctuations is that the former are cou- 
pled to long-range Coulomb interactions while the latter 
are not. This has dramatic consequences on the polariza- 
tion dependence of the observed spectra. We now elabo- 
rate upon this further. 

We limit focus to the case of non-resonant x-ray scat- 
tering since resonant processes have not yet been treated 
exactly in any correlated itinerant model. The inelas- 
tic x-ray response is given formally by a generalized 
density-density correlation function S{q,uj) — —^[1-1- 
"-('^)]x"(q,t^) with 

x(q,t^) ([p(q),/5(-q)])(^) (i) 

formed with an "effective" density operator given by 

P(q) = Yl 7a(k)ct(k + q/2)c.(k - q/2), (2) 

n{uj) is the Bose distribution function, and the double 
prime superscript denotes the imaginary part. We relate 
the inelastic light scattering vertex 7^ to the curvature of 
the energy band e(k) = t* X^jli coskj/Vd and the light 
polarizations through 

This holds in the limit of vanishing energy transfers, 
but can also be generalized in terms of Brillouin zone 
harmonics to other non-resonant cases. Here e''^ de- 
note the incident, scattered x-ray polarization vectors, 
respectively, and we have chosen units ks ~ c — h = 1 
and have set the lattice constant equal to 1. We can 
classify the scattering amplitudes by their point group 
symmetry operations. It is customary to have Aig de- 
note the symmetry of the lattice (s-wave) and Big and 
i?2g denote two of the d-wave symmetries. For any 
dimension d > 1, if we choose — (1,1,1,...) and 
= (1, —1, 1, —1, ...), then we have the Big sector, while 
— e'' = (1,1,1,...) projects out the Aig sector since the 
B2g component is identically zero for models with only 



nearest-neighbor hopping. Thus, we can cast the scat- 
tering amplitudes into a simple form: jai^O^) ~ ~^0^) 
and 7Bj^(k) ~ t* coskj(— l)^7Vd, which recovers 

the d = 2 representations of the tetragonal point group 
symmetry operations commonly used in Cu02 systems. 
We note that if we take the pure charge vertex for Aig, 

jAig = 1, then S{q,Lu) cx Im | with e the dielec- 

tric functiorki^. 

B. Weakly correlated electrons 

It is useful to review the nonresonant response for 
weakly correlated metals to determine where we expect 
to see the role of correlations emerge. For non-interacting 
electrons the effective density response is given in terms 
of a generalized Lindhard function which incorporates 
the symmetry dependence of the light scattering ampli- 
tudes 7 in the Lindhard kernel^^. In particular, in the 
limit q^ there is no low energy inelastic light scatter- 
ing (for three-dimensions) as there is no phase space to 
create electron-hole pairs and the only excitation is the 
high-energy collective plasmon. This is analogous to the 
situation of the charge susceptibility, which vanishes at a 
flnite frequency when q = because the total charge of 
the system commutes with the Hamiltonian. For flnite 
q, the particle-hole continuum gives low-energy scatter- 
ing up to a frequency ofvpq (with vp the Fermi velocity). 
When scattering off an impurity potential Vk,k' is added, 
this sharp cut-off is smeared, and scattering occurs over 
a wide range of transferred frequencies. The density re- 
sponse at small q is given by an effective density-density 
Kubo formulcii^ 

X^^(q,») = A^F "^^'_2 , (4) 

with Np the density of states at the Fermi level and 
the relaxation rate for density fluctuations having a sym- 
metry selected by light orientations labelled by L {L de- 
notes an irreducible representation of the point group of 
the crystal, such as Aig or Big for a tetragonal crystal; 
we use L = to denote the Aig sector). Expanding the 
impurity potential in terms of a complete set of basis 
functions i^lO^) yields 

l^k.k' P=E'^i(k')rL0L(k). (5) 

L 

The width and location of the peak of the response is 
given by = Tj^Iq — + Dq^, where t^^ — 2t:NfTl 
is the scattering rate that preserves charge fluctuations 
having symmetry L, and D is the diffusion constant re- 
lated to the resistivity p by an Einstein relation, D^^ = 
2e^NFP- Here we have assumed that the impurity po- 
tential is rotationally invariant and largely independent 
of momentum transfer. Thus in this case, phase space is 
already created by the impurity scattering potential for 
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anisotropic {L ^ 0) density fluctuations coupled to the 
x-rays. However, isotropic density fluctuations (L — 0) 
are governed by the continuity equation and must van- 
ish at q=0 even in the presence of an impurity potential. 
Therefore, for L ^ channels the x-ray response 

has a Lorcnzian lineshape with a peak position and width 
which grows as for momentum transfers away from the 
zone center 5 = 0, while for L — Q (^ig); there is only 
low energy scattering for finite q due to particle number 
conservation. 



III. FORMALISM WITH CORRELATIONS 

Coulomb interactions create phase-space for particle- 
hole excitations and lead to inelastic scattering even at 
q=0 for channels not having the underlying symmetry of 
the lattice. The scattering can be enhanced when the mo- 
mentum structure of the Coulomb interaction is consid- 
ered further. For example, in a material having a nested 
or slightly nested Fermi surface (FS) at some points in 
the BZ, the resulting response would be enhanced for 
polarization orientations which highlight the nested or 
nearly nested regions of the FS^°. In the case of antifer- 
romagnetic interactions which are strong for momentum 
transfers of (tt, tt) the response is appreciably modified for 
the B\g channel in two-dimensional tetragonal systems^!. 
The dispersion of these excitations can then be tracked 
as a function of light momentum transfers q just as they 
can via neutron scattering. Thus, in principle, inelastic 
x-ray scattering systematically tracks the role of correla- 
tions and the accompanying FS instabilities by exploring 
the polarization dependence and momentum transfer de- 
pendence of the resulting spectra. 

In this paper, we are interested in carrying out calcu- 
lations in which electronic correlations can be handled 
exactly in a system which can be tuned through a quan- 
tum critical point. The Falicov-Kimball model, which 
has been used to describe a variety of phenomenon in 
binary alloys, rare-earth compounds, and intermediate- 
valence materials^, contains itinerant band electrons and 
localized electrons, in which the band electrons can hop 
with amplitude^^ t* I2\fd between nearest neighbors on 
a d-dimensional hypercubic lattice and interact via a 
screened Coulomb interaction U with the localized elec- 
trons: 



t* 



2Vd 



awi 



(6) 

where cj , Ci is the spinless conduction electron creation 
(annihilation) operator at site i and = or 1 is a clas- 
sical variable for the localized electron number at site i. 
Ef and /x control the filling of the localized and conduc- 
tion electrons, respectively. We restrict consideration to 
half filling (4c,) = (m,) = 1/2. 

In this model, at half-filling, the system possesses^^ a 
non- Fermi liquid metallic ground state for U < Uc and an 



insulating state for U > Uc- The single particle density 
of states (DOS) at the Fermi level (FL) vanishes at the 
critical Uc ~ 1.5t* and the self energy develops a pole. 
As U approaches Uc from below, a pseudogap develops 
near the FL and for U > Uc the DOS evolves into lower 
and upper Hubbard bands separated at the band centers 
by U . However, the DOS is independent of temperature 
(aside from a trivial shift due to the temperature depen- 
dence of the chemical potential, if applicable) and thus it 
is not possible to determine the particle dynamics from 
the single-particle properties alone24. 

The many-body problem is solved^^ by first recogniz- 
ing that the self energy and relevant irreducible vertex 
functions are local and then mapping the local objects of 
the lattice problem onto an effective atomic problem in a 
time-dependent dynamical field A. In this procedure, we 
are interested in calculating the local Green's function, 
which is defined by 

G(r) = -Trcf%{e-P"c{T)c\Q)S{X))lZ{\) (7) 

for imaginary times t. Here Trc/ denotes the trace 
over conduction and localized electrons and Tt denotes 
the time ordering operator. The partition function is 
Z{\) = Trc/T^(exp[-/3H]S'(A)) with the evolution oper- 
ator S defined by 



S{\) = exp 



/3 rH 

dr / dT'c'{T)X{T,T')c(T') 

Jo 



(8) 



In these equations the Hamiltonian is the atomic Hamil- 
tonian, which has t* = and all time dependence is with 
respect to this atomic Hamiltonian. 

In order to determine the Green's function anywhere 
in the complex plane, we follow the iterative algorithm 
of Jarrell2&: (i) begin with the self energy E set equal 
to zero; (ii) determine the local lattice Green's function 
from the Hilbert transform 



Giz) = / depie) 



1 



z + p, - T.{z) 



(9) 



with p(e) the noninteracting DOS (a Gaussian here); (iii) 
extract the effective medium Go from G^^(z) -I- S(z) = 
Gq^(z); (iv) calculate the new Green's function from 
Giz) = (1 - wi)Go{z) + wi/[G^\z) - U]; (v) and ex- 
tract the new self energy from ^(z) — Gq^{z) — G^^{z). 
Steps (ii) through (v) are repeated until the iterations 
converge. Sometimes we need to perform weighted aver- 
ages of the iterations to attain convergence. We usually 
work with solutions that are converged to at least one 
part in lO'^. Using this algorithm, we can determine the 
Green's function and self energy either along the imag- 
inary axis, or along the real axis. These solutions are 
then employed to calculate the inelastic light scattering 
response functions. 

The inelastic light scattering is calculated by evalu- 
ating the density-density correlation function defined in 
Eq. The Bethe-Salpeter equation for the suscepti- 
bility is shown schematically in Fig. ^ Note that there 
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FIG. 1: Coupled Dyson equations for the inelastic light scat- 
tering density- density correlation functions described by the 
scattering amplitude 7a. Panel (a) depicts the Dyson equa- 
tion for the interacting correlation function, while panel (b) is 
the supplemental equation needed to solve for the correlation 
function (the difference in the two equations is the number 
of 7a factors). The symbol F stands for the local dynamical 
irreducible charge vertex. In situations where there are no 
charge vertex corrections (like Big scattering along the zone- 
diagonal) , the correlation function is simply given by the first 
(bare- bubble) diagram on the right hand side of panel (a). 



are two coupled equations, which differ by the number 
of factors of the inelastic light scattering vertex that are 
present. The solid lines denote dressed Green's functions 
in momentum space, and the symbol F denotes the local 
irreducible charge vertex. The calculation of the rele- 
vant momentum summations implied in Fig. ^ is non- 
trivial. The starting point is to determine the direction 
in momentum space that the transferred momentum q 
lies. In this contribution we consider two different di- 
rections: (i) the zone diagonal, where q — {q,q,q, ...,q) 
and (ii) a generalized "zone boundary" , where q = 
(0, q, 0, q, 0, q, 0, g) or q = {q, tt, q, tt, q, tt, ...q, tt); in all 
cases we vary < g < tt. We choose to call the wavevec- 
tor in (ii) the zone-boundary because it reduces to the 
two-dimensional zone boundary when d — 2 and it is a 
nontrivial generalization in the infinite dimensional limit. 
If, on the other hand, we examined the true infinite- 
dimensional zone boundary, where only one dimension 
has a nonzero wavevector component, then that zone 
boundary maps onto the zone center wavevector (since 
only one of the d-components is nonzero), and there is 
no dispersion. From now on we will refer to the gen- 
eralized zone-boundary direction as the zone-boundary 
direction. 

When evaluating the density-density correlation func- 
tion, we will need to evaluate momentum summations of 
the form^ 



J2 

k 3 



Si- 
2 • 



(10) 



z + fi- S(z) - e(k) z + fi- S(z) - e(k + q) 
for the Aicr sector and 



^^eos(fc, + |) 
k J 

x(-l)^- 



(11) 



1 



z + fi- S(z) - e(k) z + fi- E(z) - e(k -I- q) 

for the Big sector. In the above equations, z denotes a 
number in the complex plane. The summation can be 
evaluated by first rewriting each momentum-dependent 
Green's function as an integral of an exponential function 



1 



M - S(z) - e(k) 



(12) 

and then expanding each band-structure energy in terms 
of the summation over each component of the wavevec- 
tor. Then the integral over momentum factorizes into an 
infinite product of one-dimensional integrals. Each inte- 
gral need be expanded just to the order of 1/d, and the 
resulting terms can be exponentiated into a form that 
has a Gaussian dependence on A. The Gaussian inte- 
gral can then be evaluated directly. When we do this, 
we find that the relevant bare susceptibilities have all of 
their q-dependence summarized in the form of two scalar 
parameters 



X = hm 



d 



cos qj 



(13) 



and 



1 

X' = hm 3 y cos I 

d 

'-^ d^oo d ^ — ' 



1j 



(14) 



i=i 



In situations where the summation in Eqs. Hll|l or H12|l 
vanish, then the response function is not renormalized 
by the irreducible charge vertex, and it can be expressed 
solely in terms of the bare response function (this phe- 
nomenon was first seen for the optical conductivity^^). 
This never occurs for the Ai^ channel, but it does for 
the Sig channel when q lies on the zone diagonal. In all 
other cases, the response function is renormalized by the 
irreducible charge vertexjSS which takes the form 



1 s. 



-•m+l 



G. 



m-\-l 



(15) 



on the imaginary axis [iuJm = j7rT(2m + 1) is the 
Fermionic Matsubara frequency and ivi — 2iTTTl is the 
Bosonic Matsubara frequency]. Here = E(ja;m) is the 
local self energy on the imaginary axis and Gm = G(iix'm) 
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is the local Green's function on the imaginary axis. These 
vertex corrections are particularly crucial for the Aig 
symmetry in order to satisfy Ward identities and particle- 
number conservation. Note that the vertex corrections 
enter for the different symmetry channels away from the 
zone diagonal because at a finite momentum transfer, 
the different symmetry representations generically mix 
together. 

The strategy for determining the final forms for the re- 
sponse functions on the real axis is to first calculate the 



response functions on the imaginary axis, then replace 
Matsubara frequency summations by contour integrals 
that surround the poles of the Fermi-Dirac distribution 
function /(w) = l/[f + exp(/3u;)] with (3 = l/T. Then the 
contours are deformed to be parallel to the real axis, and 
terms that depend on the Bosonic Matsubara frequency 
as f{ijj-\-ivi) are replaced by f{uj). Finally, we analyt- 
ically continue the Bosonic Matsubara frequency to the 
real axis. This procedure was carried out in detail for the 
Raman response^*^ and will not be repeated here. 



The final formulas for the response functions are complicated integrals of functions that depend on one of six 
different bare susceptibilities. Those six bare susceptibilities are 



1 



+ /i - T,{lu) - e Vf - X 



1 f LJ + V + ^ — Yi{liJ + v) — Xe 



(16) 



Xoiuj;X,iy) = - / dep{e) 



uj + p - T,*{uj) - e yjl- X 



1 f uj -\- V + p — ll{uj ~\- v) ~ Xe 



(17) 



X' f°° 



1 p { 1^+1/ 



[uj + p~ E(w) - e]2 



^ uj-\-h'-\-}j.— Yl{uj-\-u)-~Xe 



\ 7T^ ) 



l-X^ 



Lu + fi — T,{uj) — e 



(18) 



X' 



dep{e) 



[uj + p-E*{uj)-e]^ 



1-X^ 



Lo + p — S*(u;) 



(19) 



Xo{uj;X,iy) = — dep{e) 



uj-\-i/-\-}j.~ll{uj-\-u)~Xe 



1 ijj^u^ii—TA{LO-\~i') — Xe 771 / uj^i>^^—Yj{uj^i') — Xi 



[uj + p- S(t^) - e]3 



l-X-^ 



[uj + p — S(a;) — (] 



(20) 



-1 iJ-\-iy-\-}.L—J2{uj-\-i') — Xe 771 / cj + ?y+^i— S(cj + ?y) — Xe 



(1-^2)3/2 



w -I- /Lt — S(a;) — e 



and 



Xo(t^;X,i^) = ^ ^ / dep{e) 



(21) 



1 p / cj+iy+p-S(ttj+t/)-Xe \ 

[w + Ai-S*(cj)-e]3 



1 cj+i^+/^— S(a;-f 1/) — Xe 7-1 

2 ^ TT^x? 



j+iy+fj,-'S{uj+v)-X( 



1 



[tj + ^ - S*(tj) — e] 



-1 u;+i^+//— S(cj+z/) — Xe 771 



j+v+li-'E(u+iy)-Xe 



(l-X2)3/2' 



w -|- /Li — S* (cj) — e 
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In these equations, (z) = J dep(e)/(z— e), is the Hilbert transform of the noninteracting DOS [p(e) = exp(— e^)/Y^. 
The Aig and Big responses both can be written as 



oo 



xSC'^; ^, + ^Hlg'ttg lxSl^; ^, '^)x5(^; x, v) - x?* (u^; x, ^)] 



1 + g»U-gU+'-) ^"^^'^''') ^ 



In the case of the Aig response on the zone diagonal 
q — {q,q, ...,q), we have X = cosg and X' = cos | = 
^(1 + X)/2. In the case of the Aig response on the 
zone edge we have X ^ {1 + cosg)/2 and X' = (1 + 
cosf)/2 = (1 + VX)/2 for q = {0, q,0, q, ...,0, q) and 
X (cosq - l)/2 and X' = cos |/2 = ^(l + X)/2 for 
q = {q,T:, q,T:, q, tt). In the case of the Big response on 
the zone edge, we have X = (l + cosg)/2 and X' ~ (—1 + 
cos|)/2 = (-1 + Vx)/2 for q = (0, g, 0, g, 0, g) and 

X = (cosg- l)/2 and X' ^ -cosf/2 = -^{1 + X)/2 
for q = (g, TT, g, TT, q, tt). The case of the Big response 
on the zone diagonal is much simpler, because it does not 
have any renormalizations due to the charge vertex and 
doesn't depend on X'. It becomes 

i f°° 

J — 00 

X xS(^; X, v) - [/(c^) - f{oj + iy)]xo{uJ- X,v)}. 

(23) 

At half filling, the x-ray response function has an ex- 
tra symmetry x{^.^ ^) — X*(qi for both the Aig and 
Big channels. This symmetry is straightforward, but te- 
dious to prove. First note that at half filling we have 
G(-Lj) = G*{uj) and S(-tj) = 2^ - Using these 

results, one can directly show the following six identities: 
Xo(w; X, v) = Xq{-^ - V] X, i^); xo{uj; X, v) = x5(-^ - 
^\X, v); x'Q{i^;X, v) = -Xo{-i^-v\X, v); Xoiuj;X, v) = 
-Xo{-'^-v;X^,v); xq{^^-X,v) = xl{-^ ~ v;X,v); and 
Xg(tt';X, i^) = Xq(— — v;X.,v). Now if we substitute 
bj ~uj in the integral similar to Eq. (|22|l for x(q, —i'), 
replace f{—uj) by 1 — /(w), and employ the above identi- 
ties for the different xo's, then the integral for x(q, —v) 
can be shown to be equal to x* (q> ^) plus the imaginary 
part of an integral equal to the first term in Eq. 122|) 
without the f{uj) factor. But one can show that the re- 
sulting integral is real, which proves the symmetry for 
the response function. A similar argument shows the 
bare bubble for the Big response on the zone diagonal 
also satisfies x(q, —v) = x*(q, v). 

There is a special point in momentum space, where 
the response functions become simple again. This occurs 
at the (tt, TT, tt) point, where X = —1 and X' — 0. 



In this case, Xo = ^-iid xo is proportional to xoj so the 
Bcthe-Salpeter equation factorizes, and the susceptibility 
in Eq. H22|l becomes proportional to the bare susceptibil- 
ity for any symmetry. Hence the Aig response and the 
Big response are identical at that point in the BZ. In 
fact, this result implies that polarized measurements at 
the zone corner point can immediately show the effects 
of nonlocal charge fluctuations on the inelastic x-ray re- 
sponse functions, since any deviation of the Aig signal 
from the Big signal arises from effects of nonlocal charge 
fluctuations. We feel this may be one of the cleanest 
experimental tests for the importance of nonlocal charge 
fluctuations in a correlated many-body system. 

Finally, a careful examination of Eqs. llt)H22|l shows 
that the response function depends only on X'^. Since 
the only difference for the Aig and Big responses along 
the zone edge (for — 1 < AT < 0) is a sign change in 
X', the x-ray scattering is identical for the Aig and Big 
channels along the zone edge for —1 < AT < 0. It might 
be difficult to locate the relevant path in the BZ that 
would show this behavior in two or three dimensions, 
so examining the zone corner for the effects of nonlocal 
charge fiuctuations still remains the best bet. 



IV. RESULTS 

A. Correlated metal 

With the summary of what we expect for weakly corre- 
lated metals in mind, we present the results for U = t* /2 
at different temperatures in Figs. [21 and |21 for Big 
and Aig inelastic x-ray scattering, respectively, as a func- 
tion of transferred energy for different momentum trans- 
fers throughout the BZ measured by the momentum- 
space parameter X. Panel (a) for Figs. [3 and 13 
refers to scattering along the zone diagonal X — cosg 
for the zone-diagonal wavevector q = {q,q,q, ...,q), and 
panel (b) refers to scattering along the generalized zone 
edge [here we have q = (g, 0, g, 0, g, 0) for 1 > X = 
(1 -I- cos(7)/2 > and q = {n,q,n,q, ...,Tr,q) for > 
AT = (— 1 -I- cos q) 1 2 > — 1 The curves have been shifted 
vertically for clarity. The lowest set of curves A" = 1 
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FIG. 2: Inelastic x-ray scattering response for U = 0.5t* 
in the Big channel along (a) the BZ diagonal and (b) along 
the zone edge for the half-filled Falicov-Kimball model on a 
hypercubic lattice. The solid, dotted, short-dashed and long- 
dashed curves correspond to temperatures T = 0.1, 0.25, 0.5, 
1.0, respectively. The curves have been offset for clarity. 



FIG. 3: Inelastic x-ray scattering response U — 0.5t* in 
the Aig channel along (a) the BZ diagonal and (b) along the 
zone edge. The solid, dotted, short-dashed and long-dashed 
curves correspond to temperatures T = 0.1, 0.25, 0.5, 1.0, 
respectively. 



corresponds to Raman scattering with optical photons**. 

For the Big channel (Fig. |2Jl, a well defined low en- 
ergy Fermi-like coherence peak (below U — Q.5t*) moves 
to higher energies and broadens as one moves away from 
the zone center {X = 1), as would be expected of Lan- 
dau damping via particle-hole creation at larger qiSi^ 
(recall the Falicov-Kimball model is not a Fermi liquid 
when J7 ^ 0, but can be viewed as a "dirty" Fermi liq- 
uid for small enough U). In addition, the peak sharpens 
with decreasing temperature as the channels for Landau 
damping are lost. No particular signature can be seen 
at the energy transfer of U since it falls within the line- 
shape of the Fermi-like coherence peak, and thus the role 
of electronic correlations, while present, are obscured by 
the larger Landau damping. 

In Fig. 0]we plot the position and width of the low 
energy peak for the Big channel for momentum trans- 
fers along the BZ diagonal. The peak moves to higher 
frequencies from ~ 0.25t* for small momentum trans- 
fers X < 1 and reaches a maximum ~ t* for momentum 
transfers slightly greater than (7r/2,7r/2, • • •) before soft- 
ening as the BZ corner (tt, tt, • • •) is approached. In fact, 
the peak position is comparable to U for large momentum 
transfers in all directions. The width of the peak grows 
continually with increasing q as more and more phase- 
space is created by which charge excitations may relax. 
The width initially grows like for momentum transfers 
away from the BZ center and then slows its growth rate 
farther from the zone center (recall that X = cos q so an 
initial dependence translates into a linear dependence 
on X). 

An important difference is that the Aig results have 
no low-energy spectral weight for q = as a result of 
particle-number conservation. In a model with long- 
range Coulomb interactions, the only excitation would be 



a high energy plasmon which is soft for uncharged sys- 
tems but is pushed up to higher energies via the Higgs 
mechanism by the Coulomb interaction. In our short- 
range model, a mild peak occurs on the energy scale of 
both U and the bandwidth at the zone center (we cannot 
differentiate which one dominates). The vertex correc- 
tions do not completely remove low energy scattering for 
any finite value of q, and the low energy spectral weight 
grows for increasing q either along the zone diagonal or 
zone edge. For large q, the Aig spectra have a tempera- 
ture dependence similar to the Big response, dominated 
by particle-hole excitations. In fact, the Aig and Big 
responses are identical at the {n,Tr, ...,n) point X = —1 
due to the local approximation. Any variation in the sig- 
nal at the zone corner in different symmetry channels is 
due to nonlocal many-body correlations. 

For low q however (such as X = 0.5), the tempera- 
ture dependence is nonmonotonic due to a competition 
between increased vertex corrections, which deplete spec- 
tral weight, and decreased particle- hole damping, which 
aggregates spectral weight into the Fermi-like coherence 
peak as the temperature is reduced. It is important to 
note that for an unpolarized (partially polarized) mea- 
surement, the x-ray response is a (weighted) superposi- 
tion of the Big and Aig spectra. However, the spectra 
at small q in a metal would largely have contributions 
from the Big channel due to the significant phase space 
reduction in the Aig channel. 



B. Near Critical dynamics 

Now we turn to our results for a near critical value 
oi U = 1.5t* « Uc where the density of states vanishes 
at the Fermi level and the system undergoes a metal- 
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FIG. 4: The position (circles) and the width (diamonds) 
of the low energy peak in the Big channel for momentum 
transfers along the BZ diagonal for T — O.lt* as shown in 
Fig. 13 



insulator transition (our choice for U lies just on the in- 
sulating side of the metal- insulator transition). We plot 
in Figs. O and El the results for the Big and Aig chan- 
nels, respectively, for the same temperature ranges as in 
the previous plots. The effect of electronic correlations 
is clearly visible. For both the Big and the Aig spectra, 
two peaks become discernable at small q: the low energy 
peak (similar to the one seen for smaller values of U), 
and a non-dispersive high-energy peak (at an energy of 
roughly U corresponding to transitions between the lower 
and upper Hubbard band). Indeed the results at large q 
are more similar to the small U results since the Landau 
damping pushes the low frequency peak into the high fre- 
quency peak and further smears both peaks. Again, the 
low energy peak is removed near the zone center for the 
Aig channel, but in this case the charge-transfer peak (at 
a frequency near U) remains. 

It is important to note that even though the system is 
near critical, low energy spectral weight is visible, partic- 
ularly in the Big channel. We focus now on the spectral 
weight in this region as a function of temperature, shown 
in Fig. [71 In this low frequency region one can clearly 
see for the Big channel that the low energy spectral 
weight increases with increasing temperature throughout 
the BZ. This is most clearly seen at q = 0. For the Aig 
channel the same behavior is masked by the role of ver- 
tex corrections which reduce the spectral weight for mo- 
mentum transfers near the BZ center. Nevertheless, the 
growth of intensity with increasing temperature is clearly 
seen in both channels. The growth is particularly clear at 
low frequency transfers and for increasing transfers the 
effect vanishes and crosses over at larger frequencies to a 
region where spectral weight depletes as temperature is 
increased. The point separating these regions occurs at 
a crude isosbestic point near O.bt*, where the spectra 
are roughly independent of temperature. The isosbestic 
point becomes less well-defined for momentum transfers 
away from the zone center, and therefore it is most clearly 
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FIG. 5: Inelastic x-ray scattering response U — 1.5t* in 
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observable in Raman measurements in the Big channel. 
As the temperature is increased further, the isosbestic 
behavior disappears. 

C. Insulator 

Turning to the insulating phase, our results for U — 2t* 
for the Big and Aig channel are shown in Figs. (HI and 
O respectively. Clearly two features can be resolved in 
both the Big and Aig channel: a small, dispersive low- 
energy peak for frequencies ~ t* and a large, dispersion- 
less charge-transfer peak ~ U well separated from the 
low-energy peak. Here we see more clearly the devel- 
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FIG. 7: Detail of the low energy inelastic x-ray scattering 
response U — 1.5t* along the zone diagonal for (a) the Big 
channel and (b) the Aig channel for the temperatures shown 
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opment of the transfer of spectral weight from low fre- 
quencies to higher frequencies as temperature is lowered, 
with a more clearly defined isosbestic point separating 
the low- and high-energy transfers, as shown in Fig. IIUI 
The low-energy depletion of spectral weight and con- 
comitant increase of spectral weight at high energies 
above the isosbestic point (as T is reduced to zero) was 
recently discussed in Ref. |lj for U = At* (which lies deep 
on the insulating side of the transition) where the isos- 
bestic point is more clearly observed. The low-energy 
feature in the insulating phase is determined by ther- 
mally generated double occupancies which become un- 
populated at lower temperature. The high-energy peak 
reflects the energy scale for excitations across the Mott 
gap and is relatively dispersionless due to the local nature 



FIG. 9: Inelastic x-ray scattering response U = 21* in the 
j4ig channel along (a) the BZ diagonal and (b) along the zone 
edge. The solid, dotted, short-dashed and long-dashed curves 
correspond to temperatures T = 0.1, 0.25, 0.5, 1.0, respec- 
tively. 



of the correlations. In contrast, the low-energy feature is 
a consequence of thermally generated double occupancies 
which open a low-energy band (up to energies ^ t*) able 
to scatter x-rays. The low energy peak disperses due to 
Landau damping by the thermally generated excitations, 
created in greater numbers at larger q. These excita- 
tions are frozen out for decreasing temperature and the 
low-energy intensity disappears. Only scattering across 
the Mott gap remains at an energy transfer of U . The 
charge-transfer peak for all q broadens for increasing 
temperature while the low-energy peak gains intensity 
from zero as temperature is increased, particularly in 
the Big channel. As a consequence, both Big and Aig 
possess a non-dispersive isosbestic point — a frequency at 
which the spectra are temperature independent — around 
v ~ U /2. This result agrees with our previous results for 
[/ = 4 in which the two peaks are further separated and 
the isosbestic point is more clearly observed. 

We note that even in the insulating case there is no 
spectral weight at small energy transfers for the Aig 
channel. Thus we note that regardless of the strength 
of the correlations, the Raman response (q=0) and the 
inelastic x-ray response at small q should be dominated 
by the Big response. This can be an important diagnos- 
tic tool for investigating the nature of charge dynamics 
in different regions of the BZ due to the projection of the 
Big scattering amplitude form factors compared to Aig. 



V. SUMMARY AND DISCUSSION 

In summary, we have constructed a formally exact the- 
ory for non-resonant x-ray scattering in a system which 
can be tuned across a quantum metal-insulator transi- 
tion. We focused on the polarization and momentum 
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FIG. 10: Detail of the low energy inelastic x-ray scattering 
response U = 2t* along the zone diagonal for (a) the Big 
channel and (b) the A\g channel for the temperatures shown 
in Figs. IHlandEl 

transfer dependence of the resulting spectra as a way of 
discerning the role of electron correlations. In particu- 
lar, the way in which spectral weight is transferred over 
different frequency regions as a function of temperature 
can shed important light on the strength of the electronic 
correlations, and the momentum dependence of the ob- 
served spectra can be used to determine "hot regions" 
on the FS. In general, the temperature and polarization 
dependence of the spectrum would assist in an interpreta- 
tion of observed peaks in the x-ray spectrum of correlated 
insulators for example. 

In addition we have pointed out a number of features 
which reflect the nature of the electronic correlations. 
One important finding concerns the polarization depen- 
dence of the results for momentum transfers at the BZ 
corner (tt, tt, ...tt). In a theory in which the correlations 
are purely local wc find that the response function should 
be identical at this point for both Aig and Big scatter- 
ing geometries and that any differences can be attributed 
to the importance of non-local correlations (indeed, they 
are identical for —1 < X < along the generalized zone 
boundary). In addition we have pointed out that for low 
q the full response is dominated by the Big channel which 
projects out particle-hole excitations. Thus in this limit, 
the excitations can be directly probed and tracked as a 
function of temperature. 

There is currently limited experimental data concern- 
ing the polarization and/or temperature dependence of 
the observed spectra in either correlated metals or insu- 
lators and thus many of our predictions remain open to 
experimental verification. At this stage current experi- 
ments have focused on collective excitations such as the 
plasmon^ or orbitorkii or excitations across a Mott gap 
in correlated insulators^iSiLi&. Our theory would pre- 
dict several new effects which could serve as a fingerprint 
of the role of electronic correlations in both correlated 
metals and insulators by a systematic study of the de- 



pendence on temperature and polarization orientations. 
In particular one could use x-ray scattering to elucidate 
electron dynamics near and through a quantum critical 
metal-insulator transition. 

Our theory does not address the role of resonant scat- 
tering and the connection to multiband systems. To 
capture resonance effects, detailed information is needed 
about the energy separation of the various bands as well 
as the matrix elements which couple the valence and con- 
duction bands via light scattering. For a Mott insulator 
this would include resonant transitions between the up- 
per and lower Hubbard bands as well as between exci- 
tons. Recently this has been addressed via exact diago- 
nalization studiesi^ and a spin-polaron approacbi^,, and 
its formulation for the Falicov-Kimball model is currently 
under investigation by us. A more realistic theory for res- 
onant inelastic x-ray scattering should also include reso- 
nant transitions in which the deep core hole (created by 
the incident x-ray) decays via Auger processes and must 
also include the strong perturbing effect of the core hole 
on any intermediate states (such as band states or collec- 
tive modes of the system such as plasmons or magnons) 
accessible to scattering transitionsi^. 

We have also chosen to focus on the paramagnetic 
metal to paramagnetic insulator transition. The Falicov- 
Kimball model however possesses phases containing 
charge order and phase separation. It would be extremely 
useful to examine the excitations in the ordered phases 
via Hght scattering in this model. More generally, dy- 
namical mean field theory can be used to address the 
excitations in the ordered phase of this model as well as 
the Hubbard model. 

We close with a discussion of the applicability of our 
results for the limit of large dimensions to finite dimen- 
sional systems. One important consequence of lower di- 
mensions would be that the self energy and irreducible 
vertex function will not be strictly local and the momen- 
tum dependence may crucially alter not only the for- 
malism but also the spectral evolution of the response 
as correlations are changed via doping. Particularly the 
spectra might show dispersive features which are much 
more complex than the ones we observed in these calcu- 
lations. We again note that the inelastic x-ray spectra 
for momentum transfers at the BZ corner would be very 
useful to quantify the importance of these non-local cor- 
relations. One should note, however, that the "rough- 
ness" of the Fermi surface actually simplifies as the di- 
mensionality is lowered, so the infinite-c? results already 
include many complex geometrical effects of the infinite- 
dimensional hypercubic Fermi surface. 
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